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We present a field-theoretical method to obtain consistently the equations of motion for small 
amplitude fluctuations of the order parameter directly in real time for a homogeneous, neutral BCS 
. superconductor. This method allows to study the nonequilibrium relaxation of the order parameter 

as an initial value problem. We obtain the Ward identities and the effective actions for small phase 
the amplitude fluctuations to one-loop order. Focusing on the long-wavelength, low-frequency limit 
near the critical point, we obtain the time-dependent Ginzburg-Landau effective action to one-loop 
order, which is nonlocal as a consequence of Landau damping. The nonequilibrium relaxation of 
the phase and amplitude fluctuations is studied directly in real time. The long-wavelength phase 
q , fluctuation (Bogoliubov-Anderson-Goldstone mode) is overdamped by Landau damping and the 

• relaxation time scale diverges at the critical point, revealing critical slowing down. 
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PACS numbers: 74.40.+k, 74.90.+n, 74.20.Fg 



I. INTRODUCTION 



Nonequilibrium phenomena in superconductors continue to be the focus of attention. The dynamics of Josephson 



junctions, phase slip phenomena in the dynamics of vortices and relaxation of the order parameter and supercurrents 
are few examples of the experimental effort that probe nonequilibrium aspects of superconductivity. 

Since the original work of Abrahams and Tsuneto [H there has been an ongoing effort in trying to obtain the effective 
time-dependent description of nonequilibrium phenomena from a microscopic Bardccn-Cooper-Schrieffcr (BCS) |2l 
theory Whereas the effective Ginzburg-Landau (GL) j| description in the static limit was derived by Gor'kov 
the effective time dependent Ginzburg-Landau description is still the focus of a substantial theoretical effort. There is 
• a large body of work that established the validity of a time-dependent nonlinear Schroedinger equation that describes 
{SJ \ the dynamics of the order parameter at zero temperature \% 0, 0, [l0[ |ll} |l2| . 

At finite temperature the dynamical description is complicated by the presence of Landau damping which prevents 
a local description in time because the spectral densities feature branch cuts that prevent a derivative expansion. This 
problem was originally pointed out by Abrahams and Tsuneto jl). At finite temperature Landau damping cuts are 
unavoidable and result from processes that involve scatterings of quasiparticles in the thermal bath. In derivations of 
the effective Lagrangian for dynamical phenomena from a microscopic theory the Landau damping contribution had 
often been ignored |]l3| . 

The absorptive contributions to the effective action of long wavelength phase fluctuations at finite temperature have 
been studied by Aitchison et al. jl4j for a neutral BCS superconductor. These authors studied in detail the Landau 
damping contributions to the effective action of phase fluctuations and concluded that for temperatures < T < 0.6 T c 
the effective propagator for the phase fluctuation (Bogoliubov-Anderson-Goldstone mode) can be well approximated 
' by simple quasiparticle poles at complex energy and describe damped excitations with a linear and temperature 
^ , dependent dispersion relation and narrow widths. 

An alternative approach to study nonequilibrium aspects of superconductors is based on kinetic theory. Kopnin 
p5| studied the nonequilibrium dynamics of flux flow in clean superconductors but did not address the validity of the 
time dependent Landau-Ginzburg description near the critical point. Watts- Tobin et al. []l6| studied the validity of 
the Landau-Ginzburg description near the critical point for dirty superconductors where relaxational processes are 
dominated by (elastic) collisions. 
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However, to the best of our knowledge, the description of the relaxational dynamics for the amplitude and the phase 
of the order parameter, as well as the validity of the Landau- Ginzburg description near the critical point in clean 
neutral superconductors had been the subject of several recent studies but has not been completely understood. The 
region near the critical temperature, where |Aq(T)| <C T with Aq(T) the finite-temperature gap (order parameter), 
is the region of validity of the Ginzburg-Landau theory. 

The interest on a deeper understanding of the time dependent effective action for long-wavelength phase fluctua- 
tions has been rekindled by several recent developments. Recently there has been a substantial effort to obtain the 
time dependent effective action of long- wavelength collective excitations associated with phase fluctuations in d-wave 
superconductors @ |lf, |c|. 

In particular these studies focused on the novel Carlson-Goldman modes which are Goldstone-like modes in 
charged superconductors that emerge near the critical temperature. While at T = the Anderson-Higgs mechanism 
combines the Goldstone and gauge fields into a gapped plasma mode, near the critical temperature a novel quasi- 
particle Goldstone-like excitation, the Carlson-Goldman mode, is present in charged superconductors. This mode is 
a superposition of the Bogoliubov-Anderson-Goldstone mode, present in neutral superconductors and the long-range 
gauge field which is screened at finite temperature. In Ref. |lg it was pointed out that the existence of this mode is 
associated both with screening and Landau damping of phase fluctuations. The importance of the nonequilibrium dy- 
namics of long- wavelength phase fluctuations has also been highlighted recently within the context of high-temperature 
superconductivity [ p2| . 

Furthermore, recent experiments in ultracold alkali atoms have demonstrated the trapping and cooling of fermionic 
alkalis, in particular 40 K and 6 Li jp3|. One goal of this present experimental effort is to observe a transition to a 
neutral fermi superfluid for fermi systems with an attractive interaction between atoms in two different hyperfine 
states [^4j. Recently the spectrum of low energy collective excitations in the collisionless regime has been studied, 
in particular focusing on the emergence of Goldstone or phase fluctuations in these neutral Fermi superfluids p5| . A 
proposal for the detection of the phase transition to a neutral Fermi superfluid in 6 Li alkalis relies on the spectrum 



of long-wavelength collective (Goldstone) excitations |26 



The interest on neutral BCS Fermi superfluids is interdisciplinary, from the current experimental efforts in Fermi 
alkalis and in mixtures with Bose alkalis p3| , to neutron superfluidity in nuclear matter and neutron stars. For a 
recent discussion on neutral Fermi superfluids and their interest in a wide variety of fields see Ref. Hence the 
study of the dynamics of phase fluctuations in neutral Fermi superfluids (or neutral BCS superconductors) continues 
to be of timely interest and of experimental relevance. 

The goals of this article: In this article we focus on the nonequilibrium real-time dynamics of phase and 
amplitude fluctuations in neutral BCS superconductors in the Ginzburg-Landau regime near the critical temperature. 
In particular we obtain the effective dynamical Ginzburg-Landau description of nonequilibrium relaxation of long- 
wavelength, low-frequency fluctuations of the order parameter near the critical point. 

While previous efforts, notably by Aitchison et al. |h], |T^], focused on the long-wavelength, low-frequency effective 
action well below the critical temperature for < T < 0.6 T c our goal is to study the critical region |A (T)| <C T c 
with Aq(T) the finite-temperature gap. Our study is different from previous attempts in several respects: (i) We 
implement the Schwinger-Keldysh formulation of nonequilibrium field theory pg| ] along with the recently introduced 
tadpole method [|3l| to obtain the equations of motion for small amplitude fluctuations of the order parameter in 
real time, (ii) The equations of motion obtained with these methods are retarded, lead to the Ward identities and 
allow to establish the retarded effective action at once. Furthermore, the equations of motion describe an initial value 
problem that allows a real-time study of relaxation and damping, (iii) We then focus on the Ginzburg-Landau regime 
|Ao(T)|/T -C I and establish the dynamical Ginzburg-Landau effective action to one-loop order for long- wavelength, 
low-frequency fluctuations. This effective action is retarded and nonlocal because of Landau damping, (iv) We study 
the time evolution of small phase and amplitude fluctuations from the equilibrium configuration and reveal directly 
in real time the effect of Landau damping. 

Main results: Implementing the Schwinger-Keldysh formulation of nonequilibrium field theory and the tadpole 
method we obtain the retarded equations of motion for small fluctuations, which in turn lead to Ward identities 
associated with (global) gauge invariance both in and out of equilibrium. From these equation of motion we obtain 
the retarded one-loop effective action which is nonlocal as a consequence of Landau damping. 

We then focus on the Ginzburg-Landau (Ao <SC T) region and study the real-time relaxation of small phase and 
amplitude fluctuations. While the spectral density for phase fluctuations features a peak that suggests a Goldstone-like 
dispersion relation, the relaxational dynamics is completely overdamped as a consequence of Landau damping. 

Far away from the Ginzburg-Landau regime at low temperatures, the spectral densities for both phase and amplitude 
fluctuations feature narrow quasiparticle peaks confirming previous results jf^, |f4]. In particular, the real-time 
relaxation of long-wavelength phase fluctuations is weakly underdamped by Landau damping. 

The article is organized as follows. In Sec. || we introduce the model and the linear response formulation to 
obtain the equations of motion. In Sec. Ill we introduce the Schwinger-Keldysh formulation in the Nambu-Gor'kov 
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formalism to study the nonequilibrium aspects of Bogoliubov quasiparticles. In Sec. IV we introduce the tadpole 
method, obtain the equations of motion directly in real time and cast them in terms of an initial value problem. We 
obtain explicitly the retarded self-energies to one-loop order and their spectral representations and obtain the one-loop 
retarded effective action. In Sec. [v| we obtain the Ward identities and discuss the static limit of the self-energies. 



In Sec. VI we obtain the effective time dependent Ginzburg-Landau description focusing on the Ginzburg-Landau 



regime and the long- wavelength, low- frequency limit. In this section we provide a thorough numerical analysis of the 



real-time evolution of the relaxation of phase and amplitude fluctuations. Section VII presents our conclusions and 
poses new directions. An appendix is devoted to an alternative derivation of the Bogoliubov transformation, which 
facilitates the Schwinger-Keldysh nonequilibrium formulation. 

II. PRELIMINARIES 
A. Neutral BCS model 

The BCS Hamiltonian of a neutral electron gas is given by 

H= J2 [ d 3 x^ a ( X ,t) f-ZJW( x ,t)- fl /'d 3 a;^}(x,t)V'I(x,t)^(x,t)V'T(x,t), (2.1) 
<r=t,l \ mj J 

where f/v(x, t) are the Heisenberg complex fields representing electrons of mass m and spin a, and g > is the strength 
of the s-wave pairing interaction between spin-up and spin-down electrons close to the Fermi surface. In this article, 
we set H = fee = 1. The field i/v(x, i) and its Hermitian conjugate satisfy the equal-time anticommutation relations 

{ W (x, t) , i>l (x' , t) } = <W <5 (3) (x - x') , 
{W(x, f),^ ( x \ <)} = tyt (x, t), i>l (x', t)} = 0. (2.2) 

The Hamiltonian H is invariant under the U(X) gauge transformation 

l/fe(x,t) -►e < V <r (x,t), 

^( x ,t)^e-^(x,i), (2.3) 

where 9 is a constant phase. A consequence of this U(l) gauge symmetry is conservation of the number of electrons. 
Indeed, the number operator of electrons 

N = E / <Px1>t(x,t)Mx,t) (2-4) 

<r=T4 

commutes with H and hence is a constant of motion. However, it is convenient to work in the grand-canonical 
ensemble in which the grand-canonical Hamiltonian is given by 

K = H-fiN 

= J2 I d 3 xtl>Ux,t)(-^--fj)Mx,t)-g I ^d 3 x^(x,t)V{(x,t)^(x,t)V T (x,t), (2.5) 

where the chemical potential [i is the Lagrange multiplier associated with conservation of number of electrons. The 
chemical potential /i is determined by fixing the number of electrons and in general is a function of the temperature. 
However, for the situation under study in which the temperature is much lower than the Fermi temperature, /j, can 
be approximated by its zero-temperature value, i.e., the Fermi energy. The Lagrangian (density) corresponding to K 
is given by 

- £ ^ OJ; + T- +m) Vv+S^I^'t- (2-6) 

Introducing the auxiliary complex scalar "pair field and its Hermitian conjugate defined as 

A(x,t) = gil>i(x,t)i/>i(x,t), 

At(x,t)=^|(x,i)Vi t (x,t) ) (2.7) 
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and performing the Hubbard-Stratonovich transformation |32| , the Lagrangian can be written as 

C[^^,AlA}= £ i>\ f l ^ + Yl + A^ + A^i^+^\A--A^A. (2.8) 
»=T4 V m / 9 

We note that the pair field A(x, t) is not a dynamical field as there is no corresponding kinetic term in the Lagrangian 

In the superconducting phase, we decompose the paring field into the condensate and noncondensate parts 

A(x,i) = (A(x,i))+ X (x,i), ( X (x,t))=0, (2.9) 

where (C(x, t)) = Tr[pO(x, i)]/Trpdenotes the expectation value of the Heisenberg operator 0(x, t) in the initial 
density matrix p, (A(x, t)) is the superconducting order parameter, and x( x > t) describes the noncondensate operator. 
The presence of the condensate (A(x, t)) ^ leads to spontaneous breaking of the U(l) gauge symmetry. In the 
absence of explicit symmetry breaking external sources, the condensate is homogeneous (i.e., space-time independent) 
(A(x, t)} — A , which is the situation under consideration in this article. 



B. Real-time relaxation in linear response 



The goal of this article is to obtain directly in real time the equations of motion for small amplitude perturbations 
of the homogeneous superconducting condensate in an initial value problem formulation. Our strategy to study the 
relaxation of the condensate perturbation as an initial value problem begins with preparing a superconducting state 
slightly perturbed away from equilibrium by applying an external source coupled to the pair field. Once the external 
source is switched off, the perturbed condensate must relax towards equilibrium. It is precisely this real-time evolution 
of the nonequilibrium fluctuations around the condensate the focus of this article. 



Let rj(y:,t) be an external c-number source coupled to the pair field A(x,t), then the Lagrangian given by (2.S) 
becomes 

£[V> f ,V,At,A] -> C[il>\i/>,A\A] + A^ + yf A. (2.10) 

The presence of the external source r\ will induce a (linear) response of the system in the form of an induced expectation 
value 



(A(x,t))„=A + *(x,t). (2.11) 

Here, (A(x, t))^ denotes the expectation value of the paring field A(x, t) in the presence of the external source, 
A = (A(x, t)) is the homogeneous order parameter in the absence the external source, and S(x, t) is the space- 
time dependent perturbation of the homogeneous condensate Ao induced by the external source. The linear response 
perturbation (5(x, t) vanishes when the external source 7j(x, t) vanishes at all times. This is tantamount to decomposing 
the field into the homogeneous condensate (Ao), a small amplitude perturbation induced by the external source 
[<5(x, t)], and the noncondensate part [x(x, t)] as 

A(x,t) = A + 5(x,t) + X (x,t), (x(x,t)>„ = 0. (2.12) 

In linear response theory <5(x, t) can be expressed in terms of the exact retarded Green's function of the pair fields in 
the absence of external source || Q . An experimentally relevant initial value problem formulation for the real-time 
relaxation of the condensate perturbation can be obtained by considering that the external source is adiabatically 
switched on at t — — oo and switched off at t = 0, i.e., 

r?(x,t) =77(x)e et 9(-t), e^0+. (2.13) 

The adiabatic switching-on of the external source induces a space-time dependent condensate perturbation <5(x, t), 
which is prepared adiabatically by the external source with a given value <5(x, 0) at t = determined by Tj(x). For 
t > after the external source has been switched off, the perturbed condensate will evolve in the absence of any 
external source relaxing towards equilibrium. Thus, the external source rj(x,t) is necessary for preparing an initial 
state at t = setting up an initial value problem. This method has been applied to study a wide variety of relaxation 
phenomena in different settings p3| , [34|, including the relaxation of condensate fluctuations in homogeneous Bose- 
Einstein condensates |35|. 



5 



Using the decomposition ( 2.12] ) we expand the Lagrangian density, and consistently with linear response, keep only 
the linear terms in S and 5* , which are the small amplitude perturbations from the homogeneous condensate induced 
by the external source 77. The Lagrangian C becomes (in the presence of the external source 77) 



£[^,ip,x^,x] = £ol^, ip,x\x] +£int[^,^,X^,x], 



(2-14) 



with 



/ d V \ 1 

<r=T4 \ m J g 

Ant^^.X^x] = (S* + X t )^T ~ -(** + A5)x + ?7*X + H.c., 



(2.15) 



where we have discarded the c-number (field operators independent) terms. We note that the Lagrangian C[ip' , ip, X > x] 
is obviously invariant under the gauge transformations 



/ ax £0/2 / £0 id a £0 x £0 



V^X 1 , A5, <5*, 77* 



-'^V^e-^xne-^Ace-^^e-^, 



(2.16) 



which, as will be seen below, is at the heart of the Ward identity. 

Whereas in general a gauge transformation is invoked to fix the condensate A to be real for convenience, this 
choice corresponds to fixing a particular gauge, which in turn hides the underlying gauge symmetry. In order to 
obtain the Ward identity associated with this symmetry we will keep a complex condensate Ao and analyze in detail 
the transformation laws of the various contributions to the equations of motion. 

The study of the static and dynamical properties of the BCS theory is simplified by introducing the Nambu-Gor'kov 
formulation. Let us introduce the two-component Nambu-Gor'kov fields j^, 36 



y> T (x, t) 

?A{(x, t) 



* T ( x ,i) = v|(x,t),</a(M) 





'0 r 




'0 0" 




'1 " 







, <7_ = 


1 


) CT 3 = 


-1 



and the 2x2 Pauli matrices 



in terms of which the Lagrangian can be written as 
with 



(2.17) 



(2.18) 



(2.19) 



x f >x] 



.a /V 2 . 



Si 



2m 



* - -x T x= 

5 



£ int [¥,t>, X \x} = 1' t ( T + (5 + x)* + ^X t (Ao + (5)+x t ?7 + H.c. 



(2.20) 



III. NONEQUILIBRIUM FORMULATION 
A. Generating functional 

The general framework to study of nonequilibrium phenomena is the Schwinger-Keldysh formulation j22| , which we 
briefly review here in a manner that leads immediately to a path integral formulation. 

Consider that the system is described by an initial density matrix p and a perturbation is switched on at a time ioj 
so that the total Hamiltonian for t > to, H(t), does not commute with the initial density matrix. The expectation 
value of a Heisenberg operator 0(t) = U^ 1 (t,to)0(to)U(t,to) is given by 

Tip 
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FIG. 1: The contour C in complex time plane in the Schwinger-Keldysh formulation. It consists of a forward branch running 
from t — — oo to t = +oo, a backward branch from t = +oo back to t — — oo, and an imaginary branch from t = — oo to 
t = — oo — i/3. The sources J serve to generate the nonequilibrium Green's functions of the Nambu-Gor'kov fields. 



where U(t, to) is the unitary time evolution operator in the Heisenberg picture 



U(t,to) = Texp 



dt'H(t') 



(3.2) 



with T the time-ordering symbol. If the initial density matrix p describes a state in thermal equilibrium at inverse 
temperature (3 with the unperturbed Hamiltonian H (t < to) = H, i.e., 



p = e 



U(t o -i0,to), 



then the expectation value (13.11) can be written in the form 



m)) 



TrU(t - i/3, h)U-\t^ t)OU{t, t ) 
TrU(to - i0,to) 



(3.3) 



(3.4) 



The numerator of this expression has the following interpretation: evolve in time from to up to t, insert the operator 
O, evolve back from t to the initial time to and down the imaginary axis in time from to to to — i/3. The denominator 
describes the evolution in imaginary time which is the familiar description of a thermal density matrix. We note that 
unlike the S'-matrix elements or transition amplitudes, expectation values of Heisenberg operators require evolution 
forward and backward in time (corresponding to the U and on each side of the operator O). 

The time evolution operators have a path-integral representation f32|| in terms of the Lagrangian, and the insertion 
of operators can be systematically handled by introducing sources coupled linearly to the fields, i.e., 



(3.5) 



where J and ,P are Grassmann-valued variables. The introduction of sources J, J\ j, and j* also allows a systematic 
perturbative expansion. In such an expansion, powers of operators are obtained by functional derivatives with respect 
to these sources, which are set to zero after functional differentiation. We note that the sources j, j* introduced 



in (3.5) to generate the perturbative expansion for the pair field s in t erms of functional derivatives with respect to 
these, are different from the external sources rj, rj* introduced in ( 2.1C ) to generate an initial value problem in linear 
response and to displace the condensate from equilibrium. 

Since there are three different time evolution operators, the forward, backward and imaginary, we introduce three 
different sources for each one of these time evolution operators, respectively. Taking to — > — oo, we are led to 
considering the generating functional [p9| 

Z[J + , J-, J?} = Trt7(-oo - i/3, -oo, J /3 )C/(-oo, +oo, J-)Z7(+oo, -oo; J+), (3.6) 

where U(tf, tn J) is the time evolution operator [see ( |3.2| )] in the presence of the source J and for simplicity of notation 
we have suppressed the spin index a and the complex conjugate of the sources JK The denominator in (3.4) is given 
by Trp = Z[0, 0, 0]. The generating functional Z[J + , J~ , J°\ can be written as a path integral along the contour in 
complex time plane (see Fig. [j]) 



Z[J+, J~, J 13 ] = / Vc^Vc^ exp 



i l d i xC c [^\^]J] 
c 



(3.7) 
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where Dq^T^c^ denotes the functional integration measure along the contour C and 



<• r+oo r+oo 

/ d i xL c [^ ] ,^J] = / d 4 x£[^+,^+; J+] - / d 4 x £[¥-,■>$>-■, J~] 

J C J — OO J — OO 



— OO — 2/3 



(3.8) 



with f_°° d 4 x = J d 3 x f_°° dt, etc. Because of the trace and the fermionic nature of the operators, the path integral 
along the contour C requires antiperiodic boundary conditions on the fields. The superscripts + and — refer to fields 
defined in the upper and lower branches, respectively, corresponding to forward (+) and backward (— ) time evolution, 
while the superscript (3 refers to the field defined in the vertical branch running down parallel to the imaginary axis. 
The negative sign in front of the action along the backward branch is a result of the fact that backward time evolution 
is determined by f7 ~ 1 ( +oo, — oo) with U the time evolution operator. The contour source J that enters in the contour 
Lagrangian Cc in (3/7) takes the values of the sources and J@ in the respective branches as displayed in Fig. [3. 

Functional derivatives with respect to the sources in the forward branch give time-ordered Green's functions, those 
with respect to the sources in the backward branch give the anti-time-ordered Green's functions, and those with 
respect to the sources in the imaginary branch give the usual imaginary-time (Matsubara) Green's functions. While 
the sources J + , J~, and J@ introduced to obtain the Green's functions via functional differentiation are different 
in the different branches, as they generate the time-ordered, anti-time-ordered and Matsubara Green's functions, 
respectively; the external source 77, the homogeneous condensate Ao, and the departure from equilibrium 5 are c- 
numbers and hence are the same in all branches. 

Writing the Lagrangian as a free and an interaction part as £ = £q +£; nt , the generating functional can be written 
as a power series expansion in the interaction part, which in turn can be generated by taking functional derivatives 
with respect to the sources J, by identifying 



s 



t± 



±1 



w 



5J±' 
S 



J 6JW 5J0' 
As a result, the full generating functional along the contour C can be written as 



Z[J] = exp ji J d 4 x£i n t,c 



l 5J^ % 5J 



Zo[J], 



(3.9) 



(3.10) 



where free field generating functional Zq[J] is given by and Q but with , J] replaced by £ [& , J] 



B. Green's functions 
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The equation of motion for the free Nambu-Gor'kov field 'J in the presence of the source J reads 

#(x,t) = -J(x,<). 
The solution of this equation of motion is given by 

4-j(x, f) = - J d 4 x'S{x - x', t - t')J{x', t'), 
where 5(x — x', t — t') is the Green's function matrix along the contour C and satisfies 

5(x - x', t - t') =6 c (t- t') <5 (3) (x - x'), 



(3.11) 



(3.12) 



(3.13) 



with 6c (t — t') the Dirac delta function along the contour J c dt' 5c(t — t') = 1. The Green's function S^x — x', t — t') 
has the form 



S*(x - x',< - t') = S*>(x - x',i - t')Q c (t - t') + S<{x - x',< - t')Q c (t' - t), 



(3.14) 
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where ®c (t — t') is the step function along the contour and S^(x — x',i — t') obey the homogeneous equations of 
motion. 

The antiperiodic boundary conditions on the fields in the path integral, a result of the trace over fermionic fields 
in (3.6), lead to the following boundary condition on the Green's function 



lim 5(x — x', <o — t) = — lim S(x — x' ,to — i/3 — t). 

to — > — oo to — > — oo 



Since along the contour to 



oo is the earliest time and to — i(3 is therefore the latest time, (3.15) entails 
- lim S > (x-x',t -if3-t'), Vi' 5 



lim 5 < (x-x',i -if) = 

tO — > — OO LQ — > — OO 

which is the Kubo-Martin-Schwinger (KMS) condition for equilibrium correlation functions |50| 
The free field generating functional Zq[J] is now obtained by writing 

tf(x,t) = *(x,t) + *j(x,t), 
* f (x,t) = $ f (x,<) + tf^(x,t), 

which leads to the result 



(3.15) 



(3.16) 



Z [J] =^ o [0]exp 



-i / d 4 x / d 4 x'j\x)S{x~x')J{x') 
Jc Jc 



(3.17) 



(3.18) 



where and hereafter x denotes the space-time coordinates (x, t) for simplicity of notation. The sou rce independent 
term Zq[0] will cancel between the numerator and the denominator in all expectation values in ( |3.l[ ). 

Furthermore, we are interested in computing Green's functions of finite real times which are defined for fields in 
the forward (+) and backward (— ) time branches but not in the imaginary branch. For these real-time Green's 
functions the contributions to the generating functional from one source in the imaginary branch and another source 
in either the forward or backward branch vanish by the Riemann-Lebesgue lemma |3l], |33| , since the time arguments 
are infinitely far apart along the contour. Therefore the contour integrals of the source terms and Green's functions 
in the generating functional factorize into a term in which the sources are those either in the forward and backward 
branches and another term in which both sources are in the imaginary branch jjlj, ||. The latter term (with both 
sources in the imaginary branch) cancel between numerator and denominator in expectation values and the only 
remnant of the imaginary branch is through the periodic boundary conditions along the full contour in the Green's 
function. 

Thus the generating functional for real-time Green's functions simplifies to the following expression, defined solely 
along the forward and backward time branches pl|, |33| 



Z[J±jt±] = exp 



+ 00 



d 4 x {C in t[-iS/SJ n , iS/SJ + ] - C int [iS/5J^-,-iS/5J-}} 



+ 00 



cx])<; -/ / ,/\r / d 4 x / [j n (x)S ++ (x-x / )J+(x , ) + J*-(x) 
xS—{x- x')J-{x') - J n {x)S+-(x - x')J-{x') - J^-(x)S-+{x - x')J + {x')] \, (3.1!)) 



with 

S ++ {x - x 1 ) = £>(x - x', t - t')@(t - t') + S<(x - x', t - t')e(t' - t), 
S—(x - x') = 5>(x - x', * - t')Q{t' -t) + S < (x- x', t - t')Q(t - t'), 
S- + {x - x') =S > (tc- x', t - 0. 

S + -{x-x , ) = S < (Tt-x',t-t / ), (3.20) 

where now — oo < t,t' < +oo and the superscripts +, — correspond to the sources defined on the forward (+) 
and backward (— ) time branches, respectively. An important issue that must be highlighted at this stage, is that 
derivatives with respect to sources in the forward (+) time branch correspond to insertion of operators pre-multiplying 
the density matrix p and derivatives with respect to sources in the backward (— ) branch correspond to the insertion 
of operators post-multiplying the density matrix. That this is so is a consequence of the fact that the density matrix 
evolves in time as U(t, ^o)po^ _1 (i) ^o) with U(t, to) the time evolution operator. 
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These four Green's functions are not independent because of the identity 

S ++ (x - x') + S—(x ~ x') - S+-{x - x') - S~ + (x - x') = 0. 



(3.21) 



The diagonal elements in S ++ (x — x') are the normal Green's functions, representing the propagation of single 
electrons, whereas the off-diagonal elements are the anomalous Green's functions, corresponding to the annihilation 
and creation of two electrons of opposite spins, respectively. 

The functions S^(x — x'), which are solutions of the homogeneous free field equation of motion, are simply related 
to the correlation functions of the free Nambu-Gor'kov fields VP, W. Indeed, taking variational derivatives of the free 
field generating functional Z [J] with respect to J , JT , one can show that 



S> b (x-x') = -i(* a (x)*l(x')), 
S< b (x-x')=i(*l(x')* a (x)), 



(3.22) 



where and hereafter a,b — 1,2 denote the Nambu-Gor'kov indices. The expectation values in the expressions above 
are in the non-interacting thermal density matrix which corresponds to the quadratic part of the Lagrangian (Hamil- 
tonian), i.e., the density matrix that describes free Bogoliubov quasiparticles in thermal equilibrium at inverse tem- 
perature (3. 

While the matrix elements S< b {x — x') can be obtained through the usual Bogoliubov transformation to the quasi- 
particle basis we present in the Appendix an alternative derivation of these correlation functions d irect ly from 
the spinor solutions of the homogeneous equations of motion. We find that the correlation functions in (3.22) in the 
infinite volume limit are given by 



d 3 k 
(2^)3 



S<(M-Oe lk ' 



2k-(x— X ) 



(3.23) 



where 



SZ„(k, t-t') = -i [[1 - n ¥ {E k )\ S ab (k) e -«fc(*-*'> + n F (E k ) S ah (k) e^'-*') 
S< (k, t-t')=i \n F {E k ) S ab {k) e -« s *(*-*') + [1 - n F (E k )} S ab (k) e iE *(*-*') 



with 



S(k) 



-u k v k 



-u k v k 
\v k \ 2 



S(k) = 



v k 



u k v k 

2 



u k v k u k 



(3.24) 



(3.25) 



In the above expressions k = |k|, np(ui) = \ j{eP^ + 1 ) is the Fermi-Dirac distribution function, u k and v k satisfying 
1 are the Bogoliubov coefficients [see (A4)], and E k is the energy of free Bogoliubov quasiparticles 



ml 



(bogolons) 



El 



Cl + |Ao| 



(3.26) 



where £/. = k 2 /2m — /i (see the Appendix). Using the relation 1 — 71f(w) 
condition j30] 

S > (k,t-if3-t / ) = -S<(k,t-t'). 

Hence, the correlation functions for the Nambu-Gor'kov fields "P, ^ that will enter in the nonequilibrium perturbative 
expansion are completely determined by ( p.22[ )-( ^26 ). 



e' np(w), one can easily verify the KMS 



(3.27) 



C. Feynman rules 



From the generating functional of nonequilibrium Green's functions (3.19), it is clear that the effective interaction 
Lagrangian relevant for the nonequilibrium calculations is given by 

C^ t [*t± , * ±, x t±, % ±] = £ int [*t+ , ( x t+ , x +] _ £. nt [$t- , X t- , X -] , (3.28) 

where the fields 5'^ ± , 5 ,± , x ^ ± and are defined on the forward (+) and backwards (— ) time branches respectively. 
Consequently, this generating functional leads to the following Feynman rules that define the perturbative expansion 
for calculations of nonequilibrium expectation values. 
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(b) 



FIG. 2: Feynman diagrams contributing to the retarded self-energies (a) En(i — x') and (b) £12(2: — x') to one-loop order. A 
solid line with arrow denotes the Nambu-Gor'kov field, a dashed line denotes the pair field, and a solid line with open (filled) 
circle denotes the condensate perturbation 8 (8*). 



(i) There are two sets of interaction vertices denned by Cf^ t [W^, V I ,± , x^, X*] 1 those in which the fields are in the 
forward (+) branch and those in which the fields are in the backward (— ) branch. There is a relative minus sign 
between these two types of vertices. 

(ii) For each kind of fields there are four types of Green's functions corresponding to correlations between fields 
defined in the forward or backward branch. 



The Green's f unct ions of the Nambu-Gor'kov fields M/ t are gi ven b y (|3.20 ) in terms of the correlation functions 
displayed in ( 3.22 ) that are completely determined by (|3.24j ), ( 3.25 ) and ( 3.26 ). 



The Green's functions of the pair fields x> X can be obtained in an analogous manner. However, due to the 
nonpropagating nature of the pair fields, the Green's functions are local in time and hence those of fields defined 
in different branches vanish identically. 

(iii) The combinatoric factors, trace over the Nambu-Gor'kov indices for fermion loops, etc., are the same as in the 
imaginary-time (Matsubara) formulation. 



IV. RELAXATION OF CONDENSATE PERTURBATIONS: AN INITIAL VALUE PROBLEM 



A. Equations of motion 



The equations of motion for the small amplitude superconductin g co ndensate perturbation S(x) induced by the 
external source r)(x) is obtained by implementing the tadpole method |pl| , {33| . This method begins by writing the pair 
fields A ± (x) in the forward (+) and backward (— ) time branches as 

A ± (x) = A + S(x)+x ± (x), (4.1) 

where A is the homogeneous condensate in the absence of external source, 5(x) is the perturbation of the condensate 
induced by the external source which vanishes in the absence of external source, and x ( x ) are the noncondensate 
part of the pair fields in the forward and backward time branches. The external source J](x) are c-number fields and 
hence taken the same value in both forward and backward branches. The strategy to obtain the equations of motion 
for small amp litude condensate perturbation S(x), S*(x) is to consider the interaction Lagrangian Z^nt^ , 'I', x j x] 



given in (2.2C) in perturbation theory and impose the tadpole condition 

(x ± (x))r l = {x^(x))r l =0 (4.2) 



order by order in the perturbative expansion, but, consistent with linear response, only keep contributions linear in 
S(x), 6*{x). 
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Using the Feynman rules described in the previous section, to one-loop order the tadpole condition (x + (0))ij 
(x^ + (0))rj = leads to the following expression 



6(x) 



+ / d 4 x'[Z u (x-x')S(x') + E 12 (x-x')S*(x')] +T(A ,A* )-r](x) 



d 4 x(x n (0)x + (x)) 



6*(x) 



+ / d^x'[^22{x-x')8*{x') + T, 21 {x-x')5{x')]+T\A Q ,Al)-rf{x) 



= 0, 

= 0, (4.3) 



where Y> a b(x — x') are the retarded self-energies of the pair fields and the tadpole T(Aq, Aq) denotes terms independent 
of the condensate perturbation 5, S* . The diagrams for £11(2; — x') and £12(2: — x') to one-loop order are depicted 
in Fig. ^, and those for £22(2 — x') and £21(2; — x') can be obtained in an analogous manner. Explicitly to one- loop 
order we find 

Enfc-a;') = -i[(^ + (x)a^ + {x)¥ + {x')a + ^+{x')) - (^ + (x)a^^ + (x)¥-{x')a + ^-(x'))] 
= —i tr [<j-S ++ (x — x')a+S ++ {x — x) — cr_S' H (x — x')a+S ^{x 1 — x)\ , 

£12(2 - x') = -i [(^ + (x)a^ + (x)^ + (x')a^ + (x')) - (¥ + (x)a^ + (x)¥-(x')a^-(x')}] 

= -itl■[<T-S ++ (x-x , )<J-S ++ (x , -x)-<t-S + -(x-x')<t-S- + {x' -x)] , (4.4) 

where tr denotes the trace over the Nambu-Gor'kov indices. The expressions for £22(2: — x ') an d £21(2; — x') can be 
obtained, respectively, from expressions in (4.4) through the replacement cr+ «-> <r_. 
The diagrams for the tadpole T(Aq, Aq) are depicted in Fig. ||, from which we find 

Ac 



r(A ,A3) = ^-(tf + (x)*-* + (x)) 
9 



Ao 
9 



+ itr [<7_5<(0)] 



(4.5) 



A simple calculation with the nonequilibrium Green's functions of the Nambu-Gor'kov fields obtained in the previous 
section shows that the tadpole to one-loop order is given by 



T(A ,AS) = A 



d 3 q 1 - 2n F (E q ) 
J2nf 2E q 



(4.6) 



where E q on the right-hand side of (4.6) is understood to be a function of |Ao| [see ( 3.26 )1 and, hereafter, the 
momentum integral J ^ ? 3 is restricted to the electronic states near the Fermi surface. 

It is customary to rewrite the integral over the momentum q as being over the energy £ (measured from the Fermi 
surface) at the expense of introducing the density of states jV(£), which is taken to be constant near the Fermi surface, 
and cutting the integral off at Iwd with ujd being the Debye energy, thus leading to 



T(A ,AS) = Ao 



-AA(0) 



E 



1 - 2n F (E)] 



(4.7) 



where E = -y/£ 2 + |Ao| 2 and Af(0) = mk-p/l-K 2 is the density of states at the Fermi surface. Setting the external 
source 77 = 0, the equilibrium condition for the homogeneous condensate T(Aq, Aq) = becomes 



Ao 



-AA(0) 



E 



[1 - 2n F (E)} 



0. 



(4.8) 



In equilibrium and below the critical temperature, the condensate Aq ^ 0. Therefore ( |4.8| ) leads to the finite- 
temperature BCS gap equation that determines Aq(T): 



9^(0) £° f [l-2n F (£0] = l. 

Introducing the space Fourier transforms for 6(x), rj(x) and T, a f,(x — x') as 

d 3 k 



(4.9) 



*(x,t) 



(2tt) 3 



«k(t) e 



ik-x 
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-AAAAA* + 



FIG. 3: Feynman diagrams contributing to the tadpole T(Ao, Aq) to one-loop order. A wiggly line denotes the homogeneous 
condensate Ao. 

etc., we find the equations of motion in momentum space to be given by 

/+oo 
dt' [E u (k, t - t') S k (t') + E 12 (k, t - t') S*_ k (t')] + gT(A , A*) 5^ (k) - gm (t) = 0, 
-OO 

/+oo 
dt' [£ 22 (k, t - t') 5*_ k (t') + £ 21 (k, * - t>) S k (t')] + gT* (A , A*) <5< 3 )(k) - gr,*_ k (t) = 0. (4.10) 
-oo 

The equations of motion obtained from the tadpole cond ition (x~(0))»? = (x = are the same as those given 

in the above expressions. While the equations of motion (|4.10|) are obtained to one-loop order, it is straightforward to 
conclude after a simple diagrammatic analysis that the structure of the equations of motion obtained above is general 
and valid to all orders in perturbation theory. 

From the explicit expressions for the self-energies or by taking complex conjugation of the equation of motion for 
(5k (t), one can show that 

£ 21 (k,t-t') = E* ! (-k,t-0, 

£ 22 (k,i-0 = £n(-M-0- (4-11) 

Furthermore, rotational and parity invariance imply that the self-energies are only functions of k. Upon expressing the 
nonequilibrium Green's functions in terms of the correlation functions S^, one finds immediately that the retarded 
self-energies have the following causal structure 

£ u (k,*-i') = £f 1 (fc,t-*0 e(i-i')> 

£i 2 (k, t-t') = Ef 2 (fc, t - t') Q(t ~ t'). (4.12) 
Befor e proceeding further, we note that the invariance of the Lagrangian under the global gauge transformation 



(2.16) implies that all the terms in the equations of motion for 5 (5*) must transform as S (5*) itself under the gauge 
transformation. This in turn entails that the normal self-energies En and E 22 must be invariant under the gauge 
transformation (2.16), while the anomalous ones Ei 2 and E 2 i must transform as Aq and Aq 2 , respectively. Thus, it 



is proves convenient to write 

EfiOM-t') = Y>$£(k,t-t') = Z D (k,t-t'), 
Ef 2 (fc, t-t') = E«* (fc, t-t') = (A /A*)E o (fc, t - tO, (4.13) 



where both E^ andEo are invariant under the gauge transformation ( 2.16| ). While rewriting the self-energies in this 



manner may seem a redundant exercise, the main point is to highlight and make explicit their transformation laws 
under the gauge transformation. This is an important aspect that needs to be addressed carefully in order to extract 
the Ward identities, an exact result of the underlying gauge symmetry to be explored below. 

The gauge invariant self-energies E^ and Eo can be written in terms of their spectral representation as 

/+oo 
duj [A(k, oj) sinu>(t - if) + iS{k, u) cos uj(t - if)] , 
-oo 

/+oo ^ 
dwA(k,u)sintj(t-t'). (4.14) 
-OO 

The symmetric (5) and antisymmetric (A, A) spectral functions are, respectively, even and odd functions of oj and 
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are to be given by 
S(k,oj) - 

A(k,w) - 

A(k,u>) -- 



d 3 g 
(2tt) 3 



1 1 ■ ■ //, ( E„) ■ ■ /„ (£■„)] ( -|- + |i ) - E,, ■ E,,\ i- ,Su- + EX 



- [n F (E q ) - np(Ep)] (Jf - jf ) [6(u - E q + E p ) + 5(u + E q - j, 



d 3 g 
(2^)3 



[l-nF(B g )-n F (£p)] 1 + 



EqEp 



[8(w -E q - E p ) - 6{lj + E q + E p )} 



[n F 



(E q ) - n F (E p )} (l - [S(oj -E q + E p ) - S(u + E q - E p )] J, 



1 y A |A | 



{2n) 3 E q Ep 



{[1 - n F (S 9 ) - n P (£7 p )][*(w -E q - E p ) - 5{u + E q + E p )\ 



+ [n F (E q ) - n F {E p )][5(u - E q + E p ) - 5(uj + E q - Ep)}}, 



(4.15) 



where p = |k + q| and £ q = q 2 /2m — fi. 

The terms proportional to [1 — n F (E q ) — n F (E p )] in the above expressions correspond to processes in which two 
quasiparticles are either created or destroyed, whereas the terms proportional to [n F (E q ) — n F (E p )] arise from Landau 
damping processes corresponding to scattering of quasiparticle excitations in the medium. 

In an experimental situation the dynamical evolution of the small amplitude condensate perturbation is studied 
by preparing a superconducting state slightly perturbed away from equilibrium by adiabatically coupling to some 
external source in the infinite past. Once the source is switched-off at time t = the perturbed condensate relaxes 
towards equilibrium and the relaxation dynamics is studied. As discussed above, this experimental situation can be 
realized within the real-time formulation described here by taking the spatial Fourier transform of the external source 
to be of the form 



Vk{t)=ri k e et @(-t), e^0 H 



(4.16) 



The e-term serves to switch on the source adiabatically from t = — oo so as not to disturb the system too far from 
equi librium in the process. If at t = — oo the system was in an equilibriu m st ate, then the condition of equilibrium 
(4.8) ensures that for t < there is a solution of the equations of motion ( 4.1C ) of the form 



Sk(t) = S k (0)e et for t < 0, 



(4.17) 



where ^k(O) is related to 7?k through the equations of motion for t < 0. The advantage of the adiabatic switching-on 
of the external source is that the time derivative of the solution ( 4.17 ) satisfies Sk(t < 0) — > as e — > + . 
Let us introduce auxiliary quantities H^ b (k,t — t') defined as 



x2 b (k,t-t') = d t ,n*(k,t-t'), 



(4.18) 



then upon using integration by parts, neglecting terms that vanish in the adiabatic limit e — > + , and taking Ao to 
be the equilibrium condensate and hence T(Ao, Aq) = 0, we find the equations of motion (4.10) become 

[i + g u^(k, o)]<5 k (i) + 3 nf 2 (fc, o)6*_ k (t) - g [ dt 1 [nf x (fc, t - H%tf) + nf 2 (fc, t - 0£. k (f )] = o= 

Jo 

[l+gUg(k,0W_ k (t)+gU^(k,0)5 k (t)-g [ dt'[llg(k,t - t')6*_ k (t') + H^(k,t - t')S k (t')] =0. (4.19) 

Jo 

The above coupled equations of motion for the condensate perturbation are now in the form of an initial value problem 
with initial conditions specified at t = and can be solved by Laplace transform. Introducing a two-component 
Nambu-Gor'kov spinor and its Laplace transform 



s k (t) 



-k(*) 



0k (s) 



Ms) 

*-k(s) 



(4.20) 



where s is the Laplace variable, one can rewrite the Laplace transformed equations of motion in a compact matrix 
form as 



G _1 (fc, s) <pk(s) = i[G- 1 (fc, s )-G- 1 (fc,O)]0 k (O). 
s 



(4.21) 
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In the above equation, G 1 (k, s) is the inverse Green's function (matrix) to one-loop order expressed in terms of the 
Laplace variable s 



G- X (k,s) = 



l + gJl D (k,a) <?(A /AS)£o(M' 
;(A5/Ao)So(A:,s) l+gZ D (k,-s) . 



where E d and Eo are the Laplace transforms of E^ and Eo , respectively, 



E D (k, s) 
E (fc, s) 



+oc dk 



-oo - is 



- [A(fc,fc )+S(A;,feo)] 



->() 



DC 



fcg — is 



A(fc,fco). 



The solution of ( |4.22|) reads 



Ms) = - [1 - G(fc, s) G" 1 ^, 0)] k (O), 



where 



G(fc,s) 



1 



D(k,s) 

with the denominator D(k, s) given by 



l + gi: D (k,-s) -ff(Ao/A5)E (M)' 
- 5 (A5/A )E o (fc,s) l+ 5 E D (fc,s) 



£>(*,«) = [l+. 9 E D (fc, S )][l+ 5 E D (A:,- S )] - [ 5 E (fc, 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



The real-time evolution of the condensate perturbation (f>k(t) with an initial value 0k (0) is now obtained from the 
inverse Laplace transform 



4k(*) = / £ Ms), 



(4.27) 



where the Bromwich contour B runs parallel to the imaginary axis in the complex s plane to the right of all the 
singularities (poles and cuts) of 0k (s) We note that there is no isolated pole in </>k(s) at s = since the residue 
vanishes. 



B. Self-energies in the long- wavelength, low- frequency limit 



In this article we are interested to study the relaxation of long- wavelength, low-frequency fluctuations of the pair 



field, hence our next task is to expand the self-energies as a function of fe, s up to (D(k , s ). With £ q — q /2m 
and p = I q ■ 



k| we use the following approximations 

£ P « + 



<$£ = v F kcos6, 
5£ + (SO 2 |A | 2 

9 



9 £„ 2 



£ 3 
^9 



(4.28) 



where is the angle between k and q. We keep only up to terms of 0(k 2 ,s 2 ) in the Laplace transform of the 
self-energies and obtain 



£d(M = --AA(O) 



1 - 2n F (E) 
2E 



1 + h 



i 



E (fc,s) 



E 2 

at(o) rv Ao 

Jo 
<9n F (-E) 



1 



2£ 4 

2 r 1 - 2n F (£) 



3s 2 \ 
2£ 2 J 



E 2 



2E 



' ' 4 (,£■- 



~ 12£ 2 

isE 
' 2v F k£ 

5£ 



In 



AE 2 ) 

is/v F k + £/E 
is/vpk — £/E 

2 



5^ 
£ 2 



t; 2 fc 2 £ 2 
6£ 4 



3^' M 



3s 2 
2£ 2 



s 

s 2 " ~ Te 2 

isE 



e) \ 2v F k$ 



In 



is/v F k + £/E \ 
is/v F k - i/E ) 



(4.29) 
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where E = i/£ 2 + | Ao| 2 and the dots stand for terms of higher order in the ratios fc/|Ao|, s/|Ao|. We note that to 
this order both £i)(fc, s) and T,o(k, s) are even functions of the Laplace variable s. This important feature leads to 
the decoupling of the phase and amplitude fluctuations of the pair field in the equations of motion, as will be shown 
below. 



C. One- loop effective action 



The full equations of motion in real time ( 4.1 0| ) and their Laplace transform around the equilibrium state (4.21) 
allow us to obtain at once the retarded effective action in Fourier space. The retarded Green's function kernel is 
obtained from G~ 1 (k, s) in ( |4.2l| ) by the analytic continuation in the complex s plane s — > —iui + + . In terms of the 
Fourier transform (retarded) of the two-component Nambu-Gor'kov spinor and the tadpole spinor 



T k (w) 



r(A ,A* ) 

T*(A ,A£) 



the retarded one-loop effective action to quadratic order in the fluctuations is therefore given by 
S[6, 5*} = S[0, 0] + J d 3 kduj [ ( /. t _ k (-w)G- 1 (fc,w)^ k H + 2^ t _ k (-^)T k H 

where G^ 1 (k,oj) = G^ 1 (k, s = —iuj + + ) and 5[0,0] is a function of Ao, Ag such that 

WM_ r(A0(AS)i M = T*(A 0) A S ). 



dA* 



dA 



(4.30) 



(4.31) 



(4.32) 



Obviously, variational derivatives with respect to 6, S* reproduce the retarded equations of motion. We identify (4.31) 
with the one-loop effective action quadratic in the fluctuations. 



V. WARD IDENTITY AND STATIC SELF-ENERGIES 



The Ward identities, a consequence of the underlying (global) gauge invariance, are an integral part of the program 
to establish a connection with the Ginzburg-Landau description. Furthermore the equations of motion must fulfill 
these for consistency. In this section we show how the Ward identities emerge directly from the method described 
above used to obtain the equations of motion. 

A straightforward diagrammatic analysis with the Feynman rules described above reveals that the generic structure 
of the equations of motion obtained via the tadpole method remains the same to all orders in pe rturb ation theory. 



When combined with the transformation properties of Ao, S and \ under the gauge transformation ( 2.16 ), this general 
form of the equations of motion allows to derive to all orders in perturbation theory the Ward identity for the tadpole 
T(Ao,AS). 

First, consider the case in which 5, n — the tadpole condition (x + (0))?j — leads to T(Ao,Aq) = 0, which 
is the equilibrium condition for the homogeneous condensate. For a space-time independent shift of the condensate 
A — * A +<5 induced by a space-time independent source 77, the tadpole condition now leads to T(A +<5, Aq+<5*) = —77, 
which upon expanding to linear order in 5 and S* becomes 

dT(A ,A* ) dT(A ,A* ) 
T(A ,A )+ ^ 6+ ^ S =- V . (5.1) 



We now compare (5.1) with the first equation of motion given in (4.1C), which has the same generic structure as the 
full equation of motion obtained to all orders in perturbation theory, for space-time independent 5 and 8*. Comparing 
the coefficients of 6, 5* we recognize immediately that 

dT(A ,A5) 1 

7TT — = - + En(fc = 0,^ = 0), 

oA g 

dT(A a ,A* a ) 



dAl 



= Ei2(fc = 0,w = 0), (5.2) 



which relate the self-energies at zero frequency and zero momentum to the derivatives of the tadpole with respect to 
the condensate and is valid to all orders in perturbation theory. 
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The second important ingre dien t and which stems from the equations of motion ( 4.1 0| ) is that under a global 
gauge (phase) transformation ( [2.1(j ) the tadpole T(A , Aq) transforms just as S, A and 77, i.e., T(e l6 A ,e~ l6 Aq) = 
e l6 T(Ao, Ag). Taking the gauge parameter 9 to be infinitesimal and comparing the linear terms in 9, we find to all 
orders in perturbation theory the Ward identity for the tadpole 



T(Ao,A* 



dT(A ,AS) 
dA 



Ac 



dT(A ,A* ) 



dA* 



A* 



(5.3) 



Therefore upon combining (5^2) and the Ward identity (jx^), we obtain an alternative statement of the Ward identity 
which is an exact relationship between the tadpole and the self-energies at zero frequency and momentum 



Ac 



1 



En (0,0) 



ASSi 2 (0,0)=T(A ,A5), 



(5.4) 



where £ Q ;,(0,0) stands for £ Q b(fc = 0, to = 0) for notational simplicity. Above the critical temperature, both the 
tadpole and the condensate vanish thus the above equation becomes a trivial identity. However, below the critical 
temperature A ^ and hence (5.4) leads to 



gT( Ap, AS) 
A 



= 1 



En(0,0)- -^£i 2 (0,0) 
^0 



(5.5) 



It is customary to choose the condensate to be real by redefining its phase via the gauge transformation ( 2.16 ). 
However, as argued above, for a condensate with an arbitrary phase, the anomalous self-energy £12 must be propor- 
tional to Aq since in the equation of motion it multiplies S* , which transforms under gauge transformations just as 
Aq . Thi s fact can be s een e xplicitly at the one-loop order in the expre ssions for the respective anomalous self-energy 



111 ([Of) as well as in ( [4.14] ) with the spectral functions given by ( |4.15j ). Since the product £i2<5* must transform just 
as S or Ao, therefore the phase Aq/ Aq cancels the phase of Ag in £12. In terms of the gauge invariant self-energies 
E^ and £q, the Ward identity (5.5) can be cast into an explicit gauge invariant form 



gT(A ,Ag) 
Ao 



l + 5 £ D (0,0)-g£ o (0,0). 



(5.6) 



We emphasize that the Ward identity (5.5) or (5.6) is an exact relationship valid in or out of equilibrium (corresponding 
to when the tadpole T(Aq, Aq) vanishes or not, respectively). In equilibrium the equilibrium condition T(Aq, Aq) = 
implies that the self-energies in the static limit must satisfy the relation 



l + ffEzj(0,0)-sEo(0,Q) 
The long-wavelength limit of the static self-energies is defined as 



0. 



Ed(0, 0) = lim lim £r»(fc, s), 

£o(0, 0) = lim lim £o(&, s), 
k— >o s— >o 



(5.7) 



(5.8) 



which, using ( 4.29| ), are found to one-loop order given by 

'1 - 2n F (E) 



E o (0,Q) = -jS/"(0) / d£ 
Jo 

fWD I A 12 

E o (0,0)=jV(0)y dt—jjjT 



2E 



E 2 



dn F (E) |A |- 
dE E 2 



l-2n F {E) dn F (E) 



2E 



dE 



(5.9) 



The opposite limit lim s ^o li m fc— »o £d,o(^j s) yields the first terms in the expressions (|5.9| ) but not the last terms 
proportional to dn F {E)/dE. This is a consequence of the non-analyticity of the Landau damping contributions to 
the self-energies. 

Using the expression for the one-loop tadp ole given by ((0]) it becomes clear that the identities ([T^) are fulfilled 
by the precise order of limits determined by ( |5.9| ), i.e., the long- wavelength limit of the static (s = 0) self-energies. 

We now show explicitly that the Ward identity fl5.6| ) is fulfilled to one-loop order. The tadpole to one-loop order is 
given by (4/7). From the explicit form of the one-loop self-energies in the static limit (5J3), find that 



E D (0, 0) - £ o (0, 0) = -N(0) -| [1 - 2n F (E)}, 



(5.10) 



17 



which is an exact relationship to one-loop order, obtained for arbitrary Ao. 
Upon collecting the above results, we find to one-loop order that 



3^(A ,A5) 



1 + . 9 £d(0,0)-<7£o(0,0) 



l-0Af(P)J -| [1 - 2n F {E)]. 



(5.11) 



Therefore, the Ward identity ( p.6|) is manifestly fulfilled to one-loop order. This is an important advantage of the 
tadpole method of nonequilibrium field theory: the equations of motion obtained at a given order in the loop expansion 
are causal and guaranteed to fulfill the corresponding Ward identities to that order. 

We highlight that whereas the Ward identity is independent of the limits, k,s — » 0, the individual self-energies 
have different limits because of t he n on-analyticity of the Landau da mpin g contribution. The relationship between 
the self-energies and the tadpole (5.2) is only valid in the static limit (5.9). 

The results of this section, namely the Ward identity which relates the self-energies to derivatives of the tadpole 
will play an important role in the derivation of the effective Ginzburg-Landau theory. 



VI. TIME-DEPENDENT GINZBURG-LANDAU THEORY 



We are now in position to establish contact with the Ginzburg-Landau description of the long- wavelength, low- 
frequency excitations near the critical point. The Ginzburg-Landau description in terms of a functional of the order 
parameter is valid near the critical region when \T - T c \ < T c , A(T) ~ [T C (T C - T)] 1/2 < T c . The link between the 
Ginzburg-Landau (GL) theory and the microscopic theory of superconductivity in the static limit was established by 
Gor'kov @. 

Identifying the expectation value of the pair field A = (A(x, t)) as the complex order parameter, the free energy in 
Ginzburg-Landau theory for an homogeneous order parameter (absence of gradients) is given by 

^(A,A*)=a(T)|A| 2 + ^|A| 4 , (6.1) 

where near the critical temperature a(T) cx (T — T c ). 

The linearized equation of motion of the small amplitude fluctuation of the condensate can be obtained by expanding 
the free energy around the equilibrium value of gap function; namely Ao and A*,. Writing A = Ao+5 and A* = Aq+5*, 
keeping only quadratic terms in the fluctuations: 

i * s rlK A », <dT „df „, d 2 T 5 2 S 2 T 5* 2 d 2 T 
H A, ^ = H Ao, A S ) + *_ + + 55*— + - — 2 + (6.2) 

The linearized equation of motion for both 5 and 5* can then be obtained by minimizing T with respect to 5* and 
5 respectively and are given by 

£ GL 5 + Y,f^6* + T GL (A , A*) = 0, 

Eg, L <r + E^ L <5 + [T GL (A , A*)}* = 0, (6.3) 

with 

T GL (A ,AS) = Ao [a(T) + b\A a \ 2 ], 

aT GL 



£n = [E^]* = a(T) + 2b\A 

Eg- = [E^]* = bA^ (6.4) 



dA 

L 

<VA, ! ," 



The gap equation obtained from the Ginzburg-Landau free energy (6.1) is 

T GL (A o ,AS)=0, (6.5) 



which determines the equilibrium value of the order parameter for 5 = 5* =0. When the gap equation (6.5) is 
fulfilled, the matrix of the second derivatives of the free energy ( |6.l[ ) with respect to 5, 5* has a zero eigenvalue with 
eigenvector determined by the relation 

5* =-^5. (6.6) 
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Obviously this eigenvector with zero eigenvalue corresponds to a phase fluctuation and is the Goldstone boson, or the 
Anderson-Bogoliubov-Goldstone mode, associated with the b roken global gauge symmetry. 

The similarity between the equations of motion (S.3) and (4.10) for an ho mog e neou s perturbation with vanishing 
external sources, as well as the similarities between the identities (6A) and (5.2), ( |5.4|) are now obvious and suggest 
the following identification 



9 



+ En (0,0) 
£12(0,0) 



T GL (Ao, Aq), 

^11 ! 



^12 ! 



(6.7) 



where the self-energies on the left-hand side of the relations are understood as the long-wavelength limit of the 



static self energies as discussed in Sec. J 
expression for the tadpole T(A , Aq) (| 



This similarity can be put on firmer footing by expanding the one-loop 
in terms of a power series expansion in A , Aq and keeping up to third 
order consistently with the Ginzburg-Lan dau expansion. 

Using the form of the tadpole given in (4.7) and the results obtained in Ref. [$7|, we find (h = fcj 



1) 



5 T(A ,A5) 



1 - gAf(0) In 



1.14w D \ 
T J 



A 



7. 9 AA(0)C(3) 
8^ 2 T 2 



A |A | 2 + O(A5), 



(6.8) 



where C, is the Riemann zeta function with £(3) = 1.202 and T c = l.Mwpe 1 /g-^( °). For T w T c the coef ficien t of 



the linear term can be expanded as gAf(0) (T — T c )/T c , therefore from the relations ( |5.2D and the definitions (4.13) we 
obtain 



1 + ^(0,0) = gAf(0) 
ffS o (0,0) 



T-T c 



T, 

7 5 AA(0)C(3) . 



7g-^(0)C(3) 

47T 2 T 2 



Anl 2 , 



n 2 T 2 



Anl 2 . 



(6.9) 



We now focus on the long- wavelength, low-frequency self-energies in the Ginzburg-Landau regime characterized by 
vpk,s <C I An I <C T. This regime corresponds to the description of the long- wavelength, low-frequency excitations 
near the critical point and will allow a consistent analysis of the real-time dynamics of (s mall) phase and amplitude 
fluctuations. We begin by writing the long-wavelength, low-frequency self-energies as [see ( |4.29| )1 



E D (k,s) = E D (0,0)+/(fc,s), 
£o(M = £o(0,0) + J(M, 



(6.10) 



with £^0(0,0) given by (5.E). Since the self-energies are dimensionless functions of their arguments, the expansion 
in terms of k and s must involve ratios of these variables and the typical scales in the integrals. There are two, widely 
separated scales in the Ginzburg-Landau regime, namely, T sa T c and |Ao| with |An|/T -C 1. The expressions for /, J 
feature energy denominators that would lead to infrared divergences if |A | is set to zero. These divergences reflect 
the presence of inverse powers of |An| in the expansions. 

The leading terms in the expansion can be extracted by rescaling 



£ = |A |z, £ = |A le, 



with e = Vz 2 + 1 and introducing the dimensionless variables 



|Ao| 



s 

vpk ' 



T ' |A |' 
in the integrals for I(k, s) and J(k, s). We obtain 

I(k,s) =Af{0)x [n 2 I a {x) + K 2 s 2 I b (s,x) +isl c (s,x)] , 



3.11) 



(6.12) 



(6.13) 



where 



Ia{x) 

h (s, x) 
I c (s, x) 



dz 



dz 



dz 



2T(x,z) 
3xe 3 

T(x, z) 
2xe 3 

C{x, z) 



2ze 



In 



2Z 2 



1 + ^ 



is + z/e 



- 1 



4e 2 



z 

1 + - 



C(x, z)z' 



3C(x 1 z) 
2e 4 



ise , 
—— In 

2z 



2e 6 ) 

is + z/e 
is — z/e 



z/e 



(6.14) 
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with 



T(x, z) = tanh (t^) , C(x, z) = cosh 2 



Similarly, for J(k, s) we obtain 



J(k, s) — jV(Q) a; [k 2 J a (ir) + K 2 s 2 Jb(s, x) + isJ c (s, x)] , 



(6.15) 



(6.16) 



where 



Ja{x) 

J b (s,x) 

Jc (S,X) 



1 

24 
1 

~8 
I c (s, cc) 



G?2 



, 7z 2 -3 , 2z 2 -l 

^l^^) 7 C(X,Z) 



tse is + z/e 

1 In — r 

2z is — z/e 



(6.17) 



An important consequence of the long- wavelength, low- frequency expansion of the self-energies in ([5.10 ) together 
with (6.14), Q6.17 ) is that to lowest order the self-energies are even functions of the Laplace variable s. This important 
aspect leads to the decoupling of the phase and amplitude fluctuations in the equations of motion, as can be seen as 
follows. 

If we wr ite the fluctuation of the order parameter around the space-time constant equilibrium solution of the tadpole 
equation (E13) in the form 



A(x,t) = A a + S{x,t) 

= [p + 5p{ Xl t)]e ie °+ lSe ^ t \ 



(6.18) 



where Sp(x, t) and S9(x, i), respectively, are identified as the amplitude and phase fluctuations of the order parameter 
A(x,t). It is convenient to introduce the following projection vectors 



1 
1 

V2 



1 

A£ 
"An 



at = 



V2 



1 

1 Ao 
'AS. 



Ao 

AS 



(6.19) 



in terms of these vectors and the Nambu-Gor'kov spinors (4.20), the phase and amplitude fluctuations can be written 
as 



or, equivalcntly, 



S9(x,t) = -J—pt.^t), 
v2A 



6 k (s) = V2A 



Po 



1 



\/2A 



(x,t), 



9 k(s) p + ^a 
Po 



(6.20) 



(6.21) 



In terms of the Laplace transforms 59k(s) and <5p k (s), the equations of motion (4.21) now become 

1 + g^D{k, s) - gT,o(k, s) S9w(s) = - (££>(&, s) — So(fe, s) 

J s L V 

9 



1 + 9^D(k,s) + gT, {k,s) 5p k (s) = 



Ei,(M)-Eo(M)JjWk(Q)i 
(Sc(fc, s) + Eo(fc, «)) - (Su(fc, 0) + So(fc, 0))] S Pk (0). 



(6.22) 



The decoupling of the equations of motion between phase and amplitude fluctuations is a direct consequence of the 
fact that self-energies are even functions of the Laplace variable to this, lowest order. In particular the equality 
Hp(k , s) = —s) to this or der g uarantees that the vectors p and a are eigenvectors of the matrix G^ 1 (k, s) in 

( |L2lD . The right-hand sides of fl6.22| ) are identified as the Laplace transform of the source terms Jg, J p , respectively, 
that serve the purpose of setting up the initial value problem. 
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H.025 




FIG. 4: The full x dependence of f(x) is plotted for < x < 5. The inset shows f(x) in the Ginzburg-Landau regime i< 1. 



In particular to lowest order in the long- wavelength expansion, and after the analytic continuation s — > — iuj + + , 
the retarded one-loop effective action for fluctuations around the equilibrium solution is given, up to a constant, by 



S cS [69, 6p] 



|A p 

2.9 



d 3 kdu {<W_ k (-w) [1 + .gS D (fc, w) - 9^o(k, to)} S9 k (u) + 2J e _ k (-a;)<J0 k (w)} 



i y d 3 fcdo, {<5p_ k (-w) [1 + gZ D (k, u) + gX (k, uj)} Sp k (u) + 2 J p ^ k (-cj)S Pk (u;)} , (6.23) 



where Ed o(^j w ) = ^D,o(k,s = —iuj + + ) and the mom entum integral mus t be un derstood to be restricted near 
t he F ermi surface. Expanding the self-energies as in ( 6.1C ) together with ( 6.14 ), ( 6.17 ) and using the Ward identities 
(5.6) valid in equilibrium, one can easily show that the phase fluctuations are Goldstone modes. Using the explicit 



expressions for I(k, s) and J(k, s), we find that the effective action for the phase fluctuation is identical to that obtained 
in Refs. 10, 14 to quadratic order, while the effective action for the amplitude fluctuation is a new contribution. 

The approach followed here, based directly on the equations of motion in real time also allows us to obtain the 
effective action for the amplitude fluctuations on the same footing. 



Phase fluctuations 



The initial value problem for the phase and amplitude fluctuations described by the equations of motions ( |6.22 ) can 
now be studied straightforwardly. For the phase fluctuation the inverse Laplace transform is obtained by integrating 
60k(s) along the Bromwich contour in the complex s plane parallel to the imaginary axis and to the right of all the 
singularities of the Laplace transform 



F(s,x) 



where 



S0 k (s) _ 1 

i0 k (0) sF(s,x) + f(x)'- 



f(x) = I a (x) - J a (x), 

F(s, x) = s 2 [lb (s, x) - Jb (s, 



(6.24) 



(6.25) 



In the Ginzburg-Landau regime x -C 1, we find from (|6.14|) and (|6.17|) that 



/(*) 



X 


r°° dz 


1 z z ~ 






tanh 


~12 


A, ^ 


I 2 2J 



+ 0(x 2 ) 



7C(3) 
48tt 2 " 



0(x 2 



(6.26) 



Figure |i| displays a numerical evaluation of f(x) with the full x dependence. It shows clearly that f{x) <C 1 and that 
f(x) is accurately described by (6.26) in the Ginzburg-Landau regime. 
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FIG. 5: Real and imaginary parts of F(a) 



To leading order in x, we can set x — > in the arguments when evaluating F(s,x). After the analytic continuation 
s — ► — iuj + + and introducing the dimensionless variable 

~Ti (6-27) 

VfK 

+ ,0) are, respectively, 



a 



we find in the Ginzburg-Landau regime that the real and imaginary parts of F(s — —ia 
given by (here we have set x — > 0) 
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(6.28) 



where z m - ln — yj a 2 /(l — a 2 ) and we have set x — > in the argument of F(s,x). These expressions agree with those 
found by Aitchison et al. (lj] in the limit x — > |3S| . 

The real and imaginary [-Fi(a)] parts of F(a) for < a < 1 are displayed in Fig. [|. The real part is a 

monotonically increasing function of a while the imaginary part only has contribution in the Landau damping cut 
— 1 < a < 1 corresponding to —Vpk < lu < vpk. 

Isolated real poles of the Laplace transform, describing stable quasiparticle excitations, correspond to the solutions 
of the following equations 

F K (a) - f{x) = 0, Fi(a)=0. (6.29) 

It is clear from (6.26) and Fig. [|that the first equation in ( |6.2S ) can be fulfilled for f(x) <C 1 only in a region where 
Fi(a) =/= 0. Therefore in the Ginzburg-Landau regime the long- wavelength, low-frequency Laplace transform (6.24) 
has no isolated quasiparticle poles, the only singularity is a branch cut in the imaginary axis —iv-pk < s < ivpk, a 
consequence of Landau damping. 

Having understood the analytic struc ture of the Laplace transform, we can now proceed to study the time evolution 
via the inverse Laplace transform (4.27) by closing the Bromwich contour wrapping around the cut —ivpk < s < iv-pk. 
We obtain 

/i , 
— p ph (a)cos(ar)<50 k (Q), (6.30) 



where r = v-pkt and the spectral density for the phase fluctuations is given by 

1 f(x)F 1 (a) 



P P h(c 



n[F R (a)-f(x)Y+F?(a) 



(6.31) 
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FIG. 6: Spectral density p p h(a) vs a in the Ginzburg-Landau regime. 



The spectral density p v \i(a) and the real-time evolution of 86\^{t) for several values of x < 1 are displayed in Figs. ^ 
and [t], respectively. Three important features are gleaned from these figures: 



(i) 



(ii) 



(hi) 



The spectral density features a sharp peak at a value a pC ak(^) that vanishes continuously as x — > 0. This 
peak would indicate quasiparticle excitations with dispersion relation oj = a pca .k(x)vFk. The group velocity 
v g (x) ~ 14 C(3)wf^/37t 3 vanishes at the critical point x = and increases linearly with x at least within the 
range < x < 0.2. 

While the spectral density is not of the Breit-Wigner type (Lorentzian or resonance) and hence a true width 
cannot be extracted unambiguously, it is clear from the figure that qualitatively the quasiparticle excitations 
have a "width" . This width vanishes at the critical point and increases monotonically with x <C 1 at least within 
the range consistent with a Ginzburg-Landau expansion. Therefore these quasiparticles will be Landau damped 
below T c and the relaxation time scale (the inverse of the "width") diverges at the critical point x = 0. This 
expectation will be confirmed below and suggests critical slowing down of long- wavelength phase fluctuations. 

The reason that despite the appearance of a sharp peak in the spectral density near the critical point the real- 
time dynamics does not reveal the oscillations associated with a "quasiparticle pole" is clear. As x — > both 
the damping rate and the group velocity v g (x) vanish in such a way that the time scale for damping is cither 
shorter or of the same order as the time scale for the oscillation. 



Furthermore, (|6.30| ) evaluated at t = leads to the sum rule 

1,1 da 



( — Pph(a) = 1, 

J — 1 ^ 



(6.32) 



which we have confirmed numerically for a wide range of x. 

It is clear that damping becomes more pronounced for larger x and while the peak would seem to lead to oscillations 
with period 2n / a vea k(x) there are no hints of oscillatory behavior in the real-time evolution. Phase fluctuations are 
strongly overdamped without featuring a propagating mode. Hence we conclude that near the critical point, in the 
Ginzburg-Landau regime x <C 1, Goldstone modes or phase fluctuations are severely damped despite the fact that 
the spectral density features a peak that would indicate a quasiparticle "dispersion relation" lu = a pca k(a;)fc. The 
damping becomes larger for larger x and is solely a consequence of collisionless Landau damping. 

For x <C 1, we find that the nonequilibrium relaxation of the phase fluctuation is very well approximated by an 
exponential (see t he logarithmic plot in Fig. [?]). The numerical analysis clearly indicates that 59^(t) ~ 59^(0) e~ lk ^ t , 
where from ( 6.26 ) and ( 6.2S ) the damping rate 7fc(x) is found to be given by 



lk{x) 



14C(3)w F fc 
3^3 : 



for 



(6.33) 



This result reveals critical slowing down since the damping rate vanishes at the critical point. Furthermore we also see 
that for fixed x the damping rate also vanishes in the long- wavelength limit, in agreement with the expectation that the 
relaxation time scale of Goldstone bosons should diverge in the long- wavelength limit. This is one of the novel results 
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FIG. 7: Real-time evolution of the phase fluctuation 89^^) vs r in the Ginzburg-Landau regime is plotted in linear scale (top) 
and logarithmic scale (bottom). 



of this study: The long-wavelength phase fluctuations are overdamped by Landau damping in the Ginzburg-Landau 
regime but the damping rate vanishes at the critical point indicating critical slowing down. 

In the Ginzburg-Landau regime and for long- wavelength, low-frequency fluctuations the noncquilibrium retarded 
Ginzburg-Landau effective action for phase fluctuations (Goldstone modes) is given to lowest order by (we have set 
Je = 0) 



S£[56]=Af(0) 



lAoj 
2T C 



d 3 kduj {S9^{-Lu){v F k) 2 [F(k,aj) - f(x)] 89^{uj)} 



(6.34) 



up to an additive constant, where f{x) and F(k,uj) = F(a = Lj/v-pk) are given by (6.26) and (3.28), respectively. 
The imaginary part of F(u>,k) originates in Landau damping. This long- wavelength, low- frequency effective action 
leads to the equations of motion for long-wavelength phase fluctuations in the linearized approximation valid in the 
Ginzburg-Landau regime. Thus it can be genuinely called the effective time dependent Ginzburg-Landau effective 
action. It is nonlocal in time as a consequence of Landau damping and describes real-time relaxation which is 
completely overdamped. The damping rate vanishes at the critical temperature thus signaling critical slowing down. 



Away from the Ginzburg-Landau regime: x > 1 



While we have focused in the Ginzburg-Landau regime x -C f , we now study the region x > 1 away from the domain 
of validity of the Ginzburg-Landau expansion mainly with the purpose of comparing our results to those obtained in 
Refs. rag 

In this case we must keep the full x dependence of the functions I aj b and J a ,b- The spectral density has the same 
form as in (6.31 ) but now with the replacement F(a) — > F(s — —ia + + ,ir), whose real [-Fr(q;, x)] and imaginary 
[i<i(ev, x)] parts can be found straightforwardly as in the previous case. In the low-temperature limit x 3> 1, we obtain 



/(*) 
F R (a,x) 

Fi(a,x) 



1 

12a-' 

a* 

4aj' 

7TCX 



Vl-a 2 e - x/VT ^ 



2 

^ [e x Ei(-x+) + e- x Vi(-x 



(6.35) 



where Ei(x) is the exponential integral function, x± — x (^-^==^ ± 1J and the exponentially small temperature 
corrections to f(x) and F^(a,x) have been neglected. The spectral density with the full x dependence is displayed 
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FIG. 8: Spectral density p p h(a) vs a away from the Ginzburg-Landau regime. 




— x= 1.0 
-- x = 2.0 
c = 5.0 



FIG. 9: Real-time evolution of the phase fluctuation S8^(r) vs r away from the Ginzburg-Landau regime. 



in Fig. |^ for x > 1 away from the regime of validity of the Ginzburg-Landau approximation. It shows clearly the 
emergence of a sharp quasiparticle peak, which for x 3> 1 is at a pca k = 1/V3 in agreement with the results of 
Refs. |l(], [l4]. The real-time evolution of phase fluctuations in this regime is displayed in Fig. [| It is clear from 
these figures that the sharp quasiparticle peak in the spectral density results in a real-time dynamics that is weakly 
underdamped by Landau damping. This is in contrast to the real-time dynamics in the Ginzburg-Landau regime, 
where Landau damping is so severe that the phase fluctuation is overdamped and hence there is no quasiparticle 
interpretation. 

Thus the real-time evolution displayed above for this case confirms the results of Refs. |h], [lj valid well below 
the critical temperature, that a narrow quasiparticle peak describes the dynamics of phase fluctuations. As in the 
previous case discussed above, the damping rate is found to be given by 



lk{x) ~ — - — e 



for x 3> 1, 



(6.36) 



where the factor v^k is again a consequence of the Goldstone nature of the phase fluctuation, leading to a vanishing 
damping rate in the long- wavelength limit. Because the effects of Landau damping are suppressed in the low- 
temperature limit (i» 1), the damping rate becomes smaller as the temperature decreases and hence the oscillatory 
behavior associated with the "quasiparticle pole" is evidenced. In this region a local time-dependent effective theory 
is a good approximate description of the nonequilibrium dynamics. On the contrary, the real-time dynamics in the 
Ginzburg-Landau regime (x <C 1) is overdamped, dominated by Landau damping and cannot be accurately described 
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by a local effective action in real time. Thus the Ginzburg-Landau dynamics is purely dissipative. 



B. Amplitude fluctuations 



For the amplitude fluctuations in the long- wavelength, low-frequency limit the equation of motion in terms of the 
Laplace transform requires the inverse propagator 1 + gS^fc, s) + <?So(&, a). In the case of amplitude fluctuations 
we expect a "gap" of order |Ao| in the spectrum of the quasiparticle excitations. While this pole is away from the 
region of validity of the long-wavelength, low-frequency approximation we can, nevertheless, obtain a qualitative if 
not a quantitative estimate of the dispersion relation for amplitude fluctuations. 

Isolated poles: Since the imaginary parts of £_d, So are nonzero only for —vpk < lu < vpk, for k ~ and uj ~ |Ao| 
the expected quasiparticle pole will be away from the continuum. We can find the position of this pole by looking for 
the solutions of 



Re l + gZ D (k,s)+gZ (k,s) = 0. 

J s=-iuj+a+ 



(6.37) 



This equation has solutions for values of o> given by uj amp (k), which determine the dispersion relation for the amplitude 
fluctuation. 

In the Ginzburg-Landau regime the gap of the spectrum can be estimated by setting k = in the expressions of the 
self-energies and keeping the lowest order [0(a 2 )] terms in the expansion of the self-energies. After rescaling variables 
in the integrals, we obtain 



1 + 5 S D (0, a) + 5 S o (0, s) ~ ff A/"(0) 



00 



i 



4e 2 |A |2 



(6.38) 



where use has been made of the gap equation. Keeping the lowest order in x <C 1 in the Ginzburg-Landau regime, 
the integrals can be done straightforwardly and we find 



1 + ,gS D (0, a) + g S (0, a) ~ <?AA(0) - 



1 



lAnl 2 



(6.39) 



which, upon the analytic continuation a — ► — iuj + + , suggests the gap of the spectrum in the Ginzburg-Landau 
regime to be given by 2-\/2|Ao|. 

Away from the Ginzburg-Landau regime and for arbitrary values of k the equation for the dispersion relation, 
( |6.37 ) must be solved numerically. Figure [lO] displays the dispersion relation for several values of x away from the 
Ginzburg-Landau regime. The values of the gap shown are consistent with the results obtained by Aitchison et al. 
p9| . We haste to emphasize, however, that the position of these single (quasi)particle poles are away from the regime 
of validity of our approximations and must only be taken as indicative and consistent with the findings of Ref. UA but 
not as accurate dispersion relations since higher orders in the ratio s/|A | will modify these results. 

This analysis is included here with the sole purposes of (i) emphasizing that there are single (quasi) particle poles 
away from the Landau damping continuum, consistent with the expectation of a gap in the spectrum of small amplitude 
fluctuations, (ii) establishing a comparison with the results of Ref. 
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and (iii) offering at least a qualitative, if not 
a reliable quantitative, discussion of the terms contributing to the time evolution of amplitude fluctuations. The 
full dispersion relations must be obtained by keeping all terms in the self-energies which will involve a substantial 
numerical effort, a task clearly beyond the scope of this article whose focus is on the Ginzburg-Landau regime in the 
long- wavelength, low-frequency limit. 

Landau damping cut: The imaginary part associated with Landau damping arises for k ^ and is nonvanishing 
only along the Landau damping cut —v-pk < uj < vpk. Focusing on the long-wavelength, low-frequency limit and in 
the Ginzburg-Landau regime, we find 

l + gE D (k,s)+gE {k,s) ~ 2 3 [S o (0,0) + isxAf(0)I c (s,x)] , (6.40) 

where we have used the Ward identity ( |5.5| ) and neglected the contributions from I a M and J a ,b, which are subdominant 
for vpk,s <C |Ao|. We note that in the case of phase fluctuations the terms with I c and J c cancel each other in the 
difference in the self-energies because I c = J c , however, for amplitude fluctuations these terms add up and furnish the 
dominant contribution. 

We now proceed to study the real-time evolution of the amplitude fluctuation by integrating Sp k (s) along the 
Bromwich contour in the complex a plane parallel to the imaginary axis and to the right of all the singularities of the 
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FIG. 10: Dispersion 
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FIG. 11: Real and the imaginary parts of H(a) 



$Pk( s ) _ 1 H(s,x) 



where 



We obtain 



Sp k (0) sH(s,x) + h' 
xh = S o (0,0)/7V(0), H(s,x) = isl c (s,x). 

5pk(t)= f — p amp (w,fc)cos(wt)(5 / 9 k (0), 

J-oo U 

where p amp (LJ, k) is the spectral density for the amplitude fluctuation 

Pamp(^, k) = sgn(u) Z{k) loI (k) S[u; 2 - U)\ (*)] + p£f p (a), 



■Landau regime. 



(6.41) 
(6.42) 
(6.43) 



cut / \ 
PampW 



■ Im 



h + H(s, x) 



(6.44) 



s— — ia+Q+ 



In the above expression the first term arises from quasiparticle pole w amp (fc) and the second term arises from the 
Landau damping cut — v-pk < uj < vpk. 
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FIG. 12: Temperature dependence of the residue for the quasiparticle pole Z(k) in the long-wavelength limit. 



The quasiparticle pole will contribute an undamped oscillatory component to the time evolution given by 

Spl ole (t) = Z(k) cos[u amp {k)t}. (6.45) 



The residue of the quasiparticle pole Z{k) is determined by 

- 1 da 



Z(k) = l- f —p c £ p (a). 
J-i a p 



(6.46) 



Figure [T^ shows the temperature dependence of Z(k) in the long-wavelength limit vpk <C |Ao|. It reveals clearly that 
in the Ginzburg-Landau regime [x -C 1) the spectral density p a mpftw) is dominated by the Landau damping cut. 

Whereas the damping rate of the quasiparticle vanishes in the long-wavelength, low-frequency approximation, we 
expect that higher order contributions, in particular the decay into a pair of Bogoliubov quasiparticles (bogolons) p7| , 
will lead to a width of the quasiparticle pole and hence a finite damping rate if the gap the quasiparticle spectrum is 
greater than 2|A |, as seems to be the case from the previous analysis. However, a more comprehensive study of the 
dispersion relation is needed before reaching a quantitative conclusion. Again this is beyond the regime of validity of 
the long-wavelength, low-frequency approximation stu died here. 

In the Ginzburg-Landau regime i< 1, we find [see ( |6.9| )1 



(6.47) 



87T 2 X ' 

and the real and imaginary parts of H(s = —ia + + , 0) to be given by (here we have set x = 0) 
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dz 
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a<i 7ra 2 

- — 5" ln_ ' 

8 a 



(6.48) 



respectively, where z m j n = yj a 2 /(l — a 2 ). The real and imaginary parts of H(a) are displayed in Fig. 11. The Landau 
damping contribution to the spectral density in the Ginzburg-Landau regime and the long- wavelength, low-frequency 
limit is accurately described by 



h\Hi(a)\ 



[h + H R (a)Y 



H?(a) 



(6.49) 
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FIG. 13: Spectral density pam P (a) vs a in the Ginzburg-Landau regime. 
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FIG. 14: Real-time evolution of the Landau damping contribution to the amplitude fluctuation 5pk Ut (t) in the Ginzburg-Landau 
regime is plotted in linear scale (top) and logarithmic scale (bottom). 



which, as displayed in Fig. |l3|, clearly reveals a sharp peak near a ~ 0. 

The real-time evolution of the Landau damping contribution to the amplitude fluctuation Sp^ lt (t) is displayed in 
Fig. While for x < 0.05 the nonequilibrium relaxation through Landau damping can be approximated by an 
exponential, clearly such is not the case for x > 0.05. Nevertheless we still find that the Landau damping relaxation 
time scale becomes longer as x — > 0, revealing that this contribution to the nonequilibrium relaxation seems to be 
critically slowed down near the critical point. 

We close this section, by providing the nonlocal Landau damping contribution to the retarded effective action for 
long- wavelength, low-frequency amplitude fluctuations in the Ginzburg-Landau regime (we have set J p = 0) 



S^ p [6p] = 2Af(0) 



d 3 kduj 5p_k(—oj) 



7C(3)|A | 



8tt 2 T 



H{uj,k) 



<5pk(w), 



(6.50) 



up to a additive constant, where H(u>,k) — H(a — ui/vpk) given by ( 6.48 ). Obviously, S , ^ p [<5/o] leads to the retarded 
equations of motion for the long- wavelength, low-frequency amplitude fluctuations. 
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VII. CONCLUSIONS 



In this article we have focused on the real-time nonequilibrium dynamics of small phase and amplitude fluctuations 
of the order parameter in neutral BCS superconductors in the Ginzburg-Landau regime near the critical point. We 
have implemented the Schwinger-Keldysh formulation of nonequilibrium field theory combined with the novel tadpole 
method to obtain directly in real time the retarded equations of motion for small fluctuations around equilibrium. 
These equations allow to extract the one-loop effective action for the long- wavelength, low-frequency phase and 
amplitude fluctuations in the Ginzburg-Landau regime, which is characterized by \Aq(T)\/T <C 1 with Ao(T) the 
finite-temperature gap. Furthermore, the retarded equations of motion can be cast as an initial value problem to 
study the relaxation of nonequilibrium fluctuations directly in real time. 

We studied in detail the relaxation of phase fluctuations within and away from the Ginzburg-Landau regime. Despite 
the fact that the spectral density features a sharp peak with a Goldstone-like dispersion relation in the Ginzburg- 
Landau regime the relaxation is completely overdamped as a consequence of Landau damping. This is consistent with 
a purely dissipative time-dependent Ginzburg-Landau equation. However, the effective action is nonlocal because of 
Landau damping. The relaxation is exponential in time with the damping rate 7fc(T) = 14£(3)vFfc|Ao(T)|/37r 3 T. The 
factor v-pk is a consequence of the Goldstone nature of the phase fluctuations. The relaxation of phase fluctuations 
near the critical point features critical slowing down, i.e., the relaxation time scale diverges at the critical point. 

Far away from the Ginzburg-Landau regime at low temperatures, the spectral density features sharp quasiparticlc 
peaks and the nonequilibrium relaxation is underdamped in agreement with the results of Refs. [f0[ |lj. Away from 
the critical region, the contribution from Landau damping is negligible. The long- wavelength amplitude fluctuations 
are severely Landau-damped near the critical region and the relaxation time scales also feature critical slowing down. 

While we have focused on the nonequilibrium dynamics of neutral BCS superconductors, as a next step we will 
apply these methods to the case of charged superconductors to study in detail the dynamics of the Carlson-Goldman 
modes as well as the effective action including gauge fields near the critical temperature. We postpone this study to 
a forthcoming article. 
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APPENDIX A: PLANE WAVE SOLUTIONS AND THE BOGOLIUBOV TRANSFORMATION 



In this Appendix we present an alternative derivation of the correlation functions for the fields ip a , ip^ a directly from 
the plane wave solut ions of the homogeneous equations of motion (i.e., in the absence of source) for the Nambu-Gor'kov 
field given by (3.11 ) 

'V 2 \ 

cr+A - 



d 



dt 



0-3 



2m 



/' 



The plane wave solution can be written in the form 

*(x,£) = $(fc)e- 4 ^ k ' x \ 
The two-component Nambu-Gor'kov spinor obeys 



$(fc) 



#(x,t) = 0. 



(Al) 



(A2) 





-A " 




'u k ' 






-AS 
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= LO 





(A3) 

where = k 2 /2m — ji. This is an eigenvalue equation with the eigenvalues given by to = ±25fc, where Ek = 
\J + |Ao| 2 . The normalization of the positive and negative energy spinors is chosen so that $( Q ^(fc) < i>^- ) (fc) = 5 a p, 
where a, j3 = 1,2 (not to be confused with the Nambu-Gor'kov indices) correspond to u> = ±E k , respectively. 
Introducing the Bogoliubov coefficients u k , Vk given by 

1/2 



Uk 



£,k + Ek 
2E k 



Vk = Uk 



£fc + Ek 



UkVk 



A* 

2E k : 



(A4) 
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and satisfying u\ + \vu\ 2 = 1, we find that the positive and negative energy spinors are given by 



lift 
~Vk 



Uk 



(A5) 



respectively. After accounting for the interpretation of negative energy soluti ons as antiparticles, one can therefore 
write the general plane wave solution of the homogeneous equation of motion (Ad) as 



1 



V 







— i(E k t — U.-x) 




-Vk 





Uk 



i(_B fc t-k-x) 



(A6) 



where V is the quantization volume. 

At the level of second quantization, one recognizes that ( A6) is the Bogoliubov transformation, in which the operator 
6 k , create Bogoliubov quasiparticles of momentum k (energy Ek) and obey the usual canonical anticommutation 
relations. The correlation functions of the Nambu-Gor'kov fields in the density matrix that describes free Bogoliubov 
quasiparticles in thermal equilibrium at inverse temperature (3 are therefore found to be given by (in the continuum 
limit) 



<* Q (x,iK(x',i')) = 



(*t(x',i')* (x,t)> = 



d 3 k 
(2^)3 

d 3 k 



[1 - n F (E k )} S ab (k) e -«K'-0 + n F {E k ) S ab {k) e ^<*-*') 
n F (E k )S ab (k) e -«s*(*-f ) + [1 - n F (E k )]S ab (k) e * B »<*-*') 



ik-(x-x' 



,,ik-(x— x') 



where n F {Ek) is the equilibrium distribution for Bogoliubov quasiparticles of momentum k 

1 



n F {E k ) = (blbk) = (4d k ) = 



(A7) 



and iS(fc), <S(fc) are given by ( 3.25|) . 
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